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k^ , The BABAR experiment has recently obtained some important results in the search for new 

gj ■ physics in leptonic and lepton flavor violating decays, exploiting the complete datasets col- 

lected at the T(45'), T(35') and T{2S) energies. In particular, new limits on the ratio 
r(T(15') -^ r+r-)/r(T(15') -^ l^'^l^~), on lepton flavor violating decays of the T(35') and 



0) 



^ ■ T{2S), and on r decays to three charged leptons or r ^ e//i7 are presented 

1 Introduction 
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fT) , Despite being o rigin ally devoted to study CP-violation, the BaBAR experiment (described in 

C^ I detail elsewhere ^U^ has recently obtained several important results also in precision tests of 

the standard model (SM) and in searches for new physics (NP) effects. This has been possible 
thanks to the huge data sample collected by BaBAR, mostly at an energy in the e+e~center-of- 

O ! mass (CM) frame equal to the mass of the T(4S') (corresponding to an integrated luminosity of 

426 fb"-*^), but also of the lower-mass T resonances (28 fb"-*^ at the T{3S) energy and 14 fb"-*^ at 
the T{2S)). Samples of data collected just below each T resonance (42 fb~^, 2.4 fb~^ and 1.3 

k>( I fb ~^ below the T(4S'), T{3S) and T(2S'), respectively) have been used as well. 

;-H ■ Here some recent results obtained by BaBAR are shown: a test of lepton universality in T(IS') 

- - ■ decays'^, representing the most precise measurement of the ratio r(T(lS') -^ T^T~)/r(T{lS) -^ 

fi'^fi~); and three results in the search for charged lepton flavor violating decays of the T{3S) 
and T(2S) resonances , as well as of r to e/fi^^ and to three charged leptonsl^. With each 
result BaBAR proves to be able to constrain the NP theoretical models proposed for the different 
processes, as it will be explained in the devoted sections. 

2 Test of lepton universality in T(IS') decays 

In the SM, the couplings of the gauge bosons to leptons are independent of the lepton flavor. 
Aside from small lepton-mass effects, the expression for the decay width T(IS') -^ l^l~ should 
be identical for all leptons, and given by : 

rT(i5H«-4a Q6^^^(l + 2^)yi-4^, (1) 

where a is the electromagnetic fine structure constant, Q^ is the charge of the bottom quark, 
Rn{^) is the non-relativistic radial wave function of the bound hh state evaluated at the ori- 



gin, Mx is the T(IS') mass and Mi is the lepton mass. In the SM, one expects the quantity 
i?^^(T(15)) = ^^(ig)^^+^- to be very close to one (in particular, i?^^(T(15)) ~ 0.992'St). 

In the next-to-miniinal extension of the SM^, deviations of i?r/i from the SM expectation 
may arise due to a hght CP-odd Higgs boson, A^ . Present data do not exclude the existence 
of such a boson with a mass below 10 GeV/c^. A^ may mediate the following processes . 

T{IS) -^ A% ^ /+/-7 or T{IS) -^ 775(1^)7, mO^S) ^ A^ ^ /+/-. (2) 

If the photon remained undetected, the lepton pair would be ascribed to the T(IS') and 
the proportionality of the coupling of the Higgs to the lepton mass would lead to an apparent 
violation of lepton universality. The deviation of R^^ from the expected SM value depends on 
Xd = cos^^tan/? (where 9a measures the coupling of the T(IS') to the A^, and tan/3 is the 
ratio of the vacuum expectation values of the two Higgs doublets) and on the mass difference 
between A^ and r]i,{lS). Assuming X^ = 12, r(775(lS')) = 5 MeV, and the measured M^^(i5)LLU^ 

the deviation of i?r/i(T(lS')) may be as large as ^ 4%, depending on the A^ massif. 

A measurement of this ratio has already been performed, with the result i?T-^(T(lS')) = 
1.02 ± 0.02{stat.) ± 0.05{syst.)^. 

This analysis focuses on the measurement oi Rr/^{T{1S)) in the decays T{3S) -^ T{1S)7t^7t~ 
with T{1S) -^ l^l- and 1 = fi,r oi the - 1.2 x 10^ T{3S) collected by BaBAR. Only r decays to 
a single charged particle (plus neutrinos) are considered, resulting in final states of exactly four 
detected particles for both the fi^ ji~ and t^t~ samples. 

The event selection is optimized using Monte Carlo (MC) simulated events. Different selec- 
tion criteria are used for the T(IS') -^ fi'^fi~ decays (D^) and the T(IS') -^ r^r~ decays (D^), 
because in the latter the presence of neutrinos in the final state leads to a larger contamination 
from the background (mainly non-leptonic T(IS') decays and e+e~ -^ r^r~ events). The final 
selection efficiency for the reconstructed decay chains, estimated from a sample of MC simulated 
events, are e^^ ^ 45% and Crr ^ 17% for the ji^ ii~ and the t^t~ final states, respectively. 

An extended unbinned maximum likelihood fit, applied simultaneously to the two disjoint 
datasets D^ and D^^ is used to extract Rr^ = _^i£r ^ ^m/^ ^ where Nsign (Ngigr) indicates 
the number of signal events in the Djj^ (D^) sample. For the Djj^ sample, a 2-dimensional 
probability density function (PDF) is used, based on the invariant dimuon mass M^-^^- and 
M^^^_ , the invariant mass of the system recoiling against the yr-pair, defined as: M^^^_ = 



n/ 



s + M^^ — 2 • ^/s • E^^, where ^/s is the e+e~ CM energy and E^^ indicates the vr-pair energy 
calculated in the CM frame. For the D^ sample, a 1-dimensional PDF is used, based on M^^^_. 
The functional forms of the PDFs describing the signal components are modeled from a dedicated 
sub-sample consisting approximately of one tenth of the Djj^ sample, then discarded from the 
final result in order to avoid any bias. The data collected below the T{3S) resonance are used to 
model the background shapes. The result of the simultaneous fit is i?^-^ = 1.006 ± 0.013, where 
the quoted error is statistical only. Figure [1] shows the projections of the fit results for the three 
variables. 

Several systematic errors cancel in the ratio. The main systematic uncertainties are related 
to the differences between data and simulation in the efficiency of event selection, the muon 
identification, and the trigger and background filters. There is also a systematic uncertainty on 
the signal and background yields due to the imperfect knowledge of the PDFs used in the fit. 
The total systematic uncertainty, obtained by summing in quadrature all the contributions, is 
estimated to be 2.2%. Including all the systematic corrections, the ratio i?^-^ is found to be'^: 

i?r/i(T(lS')) = 1.005 zb 0.013{stat.) ± 0.022{syst.). 

No significant deviation of the ratio i?^-^ from the SM expectation is observed. This result im- 
proves both the statistical and systematic precision with respect to the previous measurement^^^. 




9.47 9.475 
M'^. (GeV/c^) 



Figure 1: 1-D fit projections for M^+^- (top left) and for M^^""" (top rigfit) in tfie D^ sample, and for M^^""" 

(bottom) in the Dr sample. In each plot the dashed line represents the background shape, while the solid line is 

the sum of signal and background contributions to the fit, and the points are the data. 



Assuming values for X^, r{rii){lS)) and Mr^^i^is) as previously stated , the present measurement 
excludes an A^ with mass lower than 9 GeV/c^ at 90% of confidence level (CL). 



3 Searches for charged lepton flavor violation 

Lepton flavor violation (LFV) can occur via neutrino oscillation, but this has never been observed 
in charged processes because the tree-level contributions are suppressed to rates not achievable 
by the current experimental sensitivity. In many extensions of the SM, enhancements of these 
rates are possible, up to a detectable level, with expected branching fractions of (9(10~^-10~^). 
An observation of LFV in charged decays would be a clear signature of NP, and improved limits 
on the branching fractions of such processes further constrain the theoretical models proposed. 
BaBAR can search for charged LFV in several typologies of decays, both of the T resonances 
and of the r leptons. 



3.1 Search for charged LFV in narrow T decays 

This analysis searches for the charged LFV decays T{nS) -^ l^r^^ with I ^ e^fi and n = 2,3, 
using the - 1.2 x 10^ T(3S') and - 1.0 x 10^ T{2S) resonances collected by BaBAR. 

The signature of the signal events consists of exactly two oppositely charged particles: a 
primary lepton, identified as an electron or a muon, with momentum close to the beam energy, 
and a secondary charged lepton or pion from the r decay (along with other neutral particles 
not reconstructed). In order to suppress background events, if the r decays leptonically, the 
primary lepton and the r-daughter are required to have different flavors. Thus, for each value of 
n, four signal channels are defined, consisting of leptonic and hadronic r decay modes, and with 
an electron or a muon as primary lepton. The main sources of background come from r-pair 
production, as well as from Bhabha and //-pair events. 

The event selection consists of several requirements, related to the particle identification and 
to the kinematics of the r-daughter. The final selection efficiencies, estimated on samples of MC 
simulated events, vary in the range (4-6)% depending on the decay mode considered. 



An extended unbinned maximum likelihood fit is performed to the distribution of the variable 
X — \pi\/Eb^ that is, the momentum of the primary lepton [pi) normalized to the beam energy 
{Eb)- The signal distribution is expected to peak at x ^ 0.97, while the r-pair background 
X distribution is smooth and approaches zero as x -^ 0.97, and the Bhabha and /i-pair events 
have instead a peaking behavior at x ^ 1. PDFs are chosen for each of these components, 
using samples of data and of MC simulated events. The signal yield Nsig is extracted and 
found consistent with the no signal-hypothesis within ±1.8cr in all the signal channels. Since 
no statistically significant signal is observed, the 90% CL upper limit (UL) on the branching 
fraction B of each decay is determined, using a Bayesian technique, in which the prior likelihood 
is uniform in B and assumes that > 0. 

In the UL calculation, the systematic uncertainties affecting the measurement are also taken 
into account. The dominant contribution to the systematic uncertainty comes from an imperfect 
knowledge of the PDF shapes. The resulting ULs are summarized in Table[T]and are of (9(10~^), 
representing the first constraints on B{T{nS) -^ e^r^), while improving the sensitivity with 



respect to the previous ULs-^^ on B(T{nS) 



/.±r^) 



Table 1: 90% CL ULs on the branching fractions B for signal decays T{nS) 



Mode 


UL (lO"'^) 


B{T{2S) -^ e±r^) 
BiT{2S) -^ n^T^) 
B{T{3S) -^ e^T^) 
B{T{3S) -^ n^T^) 


<3.2 
<3.3 
<4.2 
<3.1 



± 



1 



3.2 Search for charged LFV in the decays r^ -^ e^7 and r^ ^ ji 

Another environment for LFV processes is r decay. In particular, r^ -^ /^7 (where I — e^fj) is 
a favored decay mode in several NP scenarios, with predicted branching fractions close to the 
current experimental limits. 

Besides being a S-factory, BaBAR has been usefully employed as a r-factory as well, since 
the cross sections for the production of r-pairs and fifi-pairs are comparable. This analysis uses 
the complete BABAR dataset, which corresponds to ~ 960 x 10^ r decays. 

The reconstructed events e^e~ -^ t^t~ show a clear topology, being well divided in two 
hemispheres: the signal side, containing the /^7-pair, required to have mass and energy compat- 
ible with the T mass and the beam energy, respectively; and the tag side, which is expected to 
contain a SM r decay, reconstructed in events where the r lepton goes to one or three charged 
tracks (with undetected neutrinos). The signal side is further required to contain only one 7 
with energy greater than 1 GeV and one track identified as an electron or a muon, separated by 
an angle determined by the kinematic of the process, since 7 and / are emitted back-to-back in 
the T rest frame. 

The main sources of background come from irreducible r-pair events, /^/~7 events and 
hadronic r decays with mis-identification of the charged tt. 

Signal decays are identified by two kinematical variables: the energy difference /S.E — EP^ — 
\/~s/2 and the beam-energy constrained r mass {tuec)' The distributions of events in niEC versus 
AE are shown in Figure [2l In this plane, a region is defined as rriEC ^ [1.55,2.05] GeV/c^ and 
AE G [—0.14,0.14] GeV, and used to extract from fits the expected fractions of background 
events. The number of events in the 2a signal ellipses (0 events for r^ -^ 6^7 and 2 events 
for T^ -^ M^t) are found to be compatible with the background expectation, without evidence 
for a signal. After the estimate of the systematic uncertainties, which are mainly due to the 



efficiencies of tracking, particle identification, trigger and background filters, frequentist yu/o 
CL ULs on the branching fractions B of the signal processes are calculated using the POLE 
program^^. The results are^^ 

B{t^ -^ e^7) < 3.3 X 10"^ and B{t^ -^ ^^7) < 4.4 x 10"^ 

at 90% CL, representing the most stringent limits on LFV in these decays. 



3.3 Limits on LFV in r decays to three charged leptons 



A further search for LFV processes is performed by BaBAR in the neutrinoless decay r~ -^ 
^i^t^s ^ where U — e^fi with z = 1, 2, 3, and charge-conjugate decay modes are implied throughout 
the section. All six lepton combinations consistent with charge conservation are considered. The 
data collected at the energy of the T(4S') are used, corresponding to ^ 430 x 10^ r-pairs. 

As explained for the previous analysis, e+e~ -^ r^r~ events can be clearly separated in a 
signal and a tag hemisphere. In this case, the signal side is required to contain three charged 
particles, identified as electrons or muons, according to one of the allowed combinations. In 
addition the three-lepton system must have mass and energy compatible with the r mass and 
the beam energy, respectively. The tag r lepton has instead to decay to one charged track, with 
undetected neutrinos. Therefore, the reconstructed final states consist of exactly four charged 
tracks with net charge equal to zero. Particle identification requirements and further selection 
criteria are applied in order to reject background events, which are mainly due to qq^ Bhabha and 
/x-pair events, as well as to SM r decays. The selection efficiencies, estimated on MC simulated 
samples, vary in the range (6-13)% depending on the channel considered. 

The signal extraction is performed using the distributions AE and AMec — t^ec — ^^r? 
where AE and ttiec have been defined in the previous section, and rur is the r lepton mass. 
The distributions of the events in this plane are shown for data in Figure [3l The expected 
background rates for each decay mode are determined by fitting data in a region defined as 
AMec ^ [-0.6,0.4] GeV/c^ and AE G [-0.7,0.4] GeV. In every channel no signal candidates 
are found in the signal region and 90% CL ULs are placed on the branching fractions B, using 
the technique of Cousins and Highland^^ following the implementation of Barlow . 

The systematic uncertainties affecting the measurement are mainly due to particle identifi- 
cation efficiency; minor contributions come from tracking efficiency and errors in the background 
estimation. 



After including all the uncertainties, the 90% CL ULs on B(r 
summarized in Table O These values supersede the previous BaBAR results^ 
with the latest limits placed by Belle*^*^. 




are calculated^, as 
and are compatible 



Table 2: 90% CL ULs on the branching fractions B for signal decays r 



Mode 


UL (10-«) 


B{t~ — )■ e~e+e~) 


<2.9 


B{t~ -> /i~e+e~) 


<2.2 


B{t~ -^ e'/i+e") 


< 1.8 


B{t~ -> ii~ ii^e~) 


<3.2 


B{t~ -^ n~e^ii~) 


<2.6 


B{t~ — )■ fj,~iJ,'^fj,~) 


<3.3 







T ^e e^e 



T" -> p," e"" e" 



iOZ 



T -^ e ^1+ p. 



T" — > fj.^ e" e" 



4 Pi 



T" — > e^ ^" ^" 



i;"->^"H^p." 



0.5 

A E (GeV) 

Figure 2: Data events (dots) for r ^ ej (top) and r -^ 

/i7 (bottom) decays. The 2a signal ellipses are shown, 

as well as the dark and light shadings representing the 

50% and 90% signal contours, respectively. 
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Figure 3: Data events (dots) for the six r decay channels 

after selection is applied. Signal regions are identified 

by the solid lines. The dark and light shadings represent 

the 50% and 90% signal contours, respectively. 
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